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Nucleosynthesis and the variation of fundamental couplings 
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We determine the influence of a variation of the fundamental "constants" on the predicted helium 
abundance in Big Bang Nucleosynthesis. The analytic estimate is performed in two parts: the first 
step determines the dependence of the helium abundance on the nuclear physics parameters, while 
the second step relates those parameters to the fundamental couplings of particle physics. This 
procedure can incorporate in a flexible way the time variation of several couplings within a grand 
unified theory while keeping the nuclear physics computation separate from any model-dependent 
assumptions. 

PACS numbers: 98.80.-k, 98.80.Ft 



I. INTRODUCTION 

The possible time variation of the fundamental cou- 
plings touches a basic cornerstone of our understanding 
of particle physics It is a characteristic feature 
for cosmological models of quintessence - combining the 
dependence of couplings on the value of the cosmological 
scalar field with the time variation of this "cosmon" field 
induces a time variation of the couplings with unfortu- 
nately unknown strength [E EL • 

Recent observations of quasar absorption spectra by 
Webb et. al. [E IE HI have suggested that the electro- 
magnetic fine structure "constant" might vary over cos- 
mological timescales, Aa em /a em = —0.54(12) x 10~ 5 for 
z » 2. However, other groups did exclude such a varia- 
tion with high statistical significance [E 0, ITTl . Also 
systematic effects, such as the evolution of isotope ratios 
could have an impact on these measurements. 
While the reality of a variation of a em in QSO absorp- 
tion lines is still in dispute we need to gain an overview 
of other possible effects of a variation of the fundamental 
couplings on cosmological observations. 

A central element of modern cosmology is Big Bang 
Nucleosynthesis (BBN). Actually, the bounds on the vari- 
ation of a em at z ~ 3 do not say much about the possi- 
ble size of a variation Aa em /a em at the time of BBN, 
around z ~ 10 10 . Furthermore, a major issue is the 
complex interplay of the variation of several couplings 
on the outcome of the element synthesis. While there 
are a number of recent investigations into the bounds of 
a variation of ote.m or other single parameters for BBN 
IE EE EE E3, EE IS , we will follow a different ap- 
proach which determines element abundances in a model 
independent way For a review of current limits on fun- 
damental couplings see |2l| . 

To determine light element abundances in the ab- 
sence of time varying couplings, one needs to know the 
particle masses and the reaction rates of the relevant 
nuclear processes (from laboratory experiments). The 
time evolution of the scale factor a(t) only depends on 
the number of relativistic particles or, more precisely, 
on the ratio p/T A for energy density and temperature 



during BBN. Excellent numerical codes for BBN abun- 
dance prediction exist |22l 1 2 .31 1241 125| . An essential pa- 
rameter for these computations is the baryon to pho- 
ton ratio r\. Taking the value from WMAP measure- 
ments ,26i], r\ — 6.14 ± 0.25 x 10~ 10 , yields a predic- 
tion of the helium abundance Y Hc = 0.2484±g;™^ [27|. 
The observational determination varies among different 
groups. Izotov and Thuan 28] quote two different val- 
ues for two different equivalent width samples of spec- 
tra. The one is Y Hc = 0.2421 ± 0.0021 and the other 
F He = 0.2444 ± 0.0020. If we were to calculate rj from 
those samples we would get n 



3.4±°;£ x 10- 10 for the 



first and rj = 4.0_q 5 x 10 for the second quote. An- 
other estimate for Yhc was obtained by Fields and Olive 
jlU F He = 0.238 ± 0.002 ± 0.005, where the second is the 
systematic error (not quoted by Izotov and Thuan). 

These results show some tension between theory and 
observation. This discrepancy is likely due to systematic 
errors which are not fully understood. For instance the 
assumptions made about the extragalactic HII regions 
differ among several groups. 

Increasing the number of light species which are effec- 
tive at BBN (e.g. more neutrinos) would enhance Yke 
and only worsen the discrepancy. (This also holds for 
the possible presence of early dark energy HEE1-) If a 
mechanism for decreasing Yee has to be found the time 
variation of fundamental couplings seems to be a partic- 
ularly plausible candidate |5| . Effects of the variation of 
the weak and strong scales or some dimensionless cou- 
pling on BBN have been discussed long ago |E HE HE| • 
One may therefore try to estimate the allowed variation 
of couplings at a very early time in cosmology. 

Confidence limits on the variation of couplings or pa- 
rameters in the framework of BBN always assume an 
underlying model. However, the confidence regions de- 
termined from a model where only a em varies are mean- 
ingless if one wants to employ a model where other cou- 
plings, such as the weak scale, are allowed to change. 
In a Grand Unified Theory (GUT) framework, not only 
does the electromagnetic interaction vary, but also weak 
and strong interactions. The details of how these are 
connected depend on the specific GUT and the varia- 
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tion of the unified couplings and mass scales of sponta- 
neous symmetry breaking. The present BBN limits on 
time varying couplings are difficult to compare due to 
this strong implicit model dependence. It is therefore 
essential to formulate the BBN estimates in a way that 
is as model independent as possible. This should facil- 
itate the comparison between different assumptions on 
the time variation of fundamental couplings. 

Our concept to solve this problem is as follows. First 
we describe our general assumption, namely that the de- 
viations are small, and explain how we linearize the prob- 
lem (Section [n}. Rather than using a numerical compu- 
tation we will use an analytic approximation to determine 
the variation of Yne ~~ this will help us to better under- 
stand how Yfjo depends on fundamental couplings. We 
emphasize that relative errors for the relative variation 
below 10 % are acceptable in contrast to the much higher 
required precision for the total abundance. In the first 
step of our analysis we estimate (Section llllll how the re- 
sults of a BBN calculation depend on seven parameters 
Xi characterizing nuclear physics (Eq. 10}). 

In a second step (Section llV|) we determine the depen- 
dence of the parameters Xi on the relevant dimcnsionlcss 
"fundamental" particle physics parameters Gk (Eq. J5J), 
which characterize the standard model of electroweak, 
strong and gravitational interactions. The connection be- 
tween both is formulated in form of a "transfer matrix" 
fik ■ The advantage of this separation is the possibility to 
compute fik without invoking BBN whereas the first step 
does not use any assumption about the particle physics - 
nuclear physics connection. The two issues can therefore 
be dealt with independently. Any improvement on the 
estimate of the dependencies in the first step can be prop- 
agated to the fundamental couplings without repeating 
the whole calculation. In the same spirit one may discuss 
a third step which relates the standard model parameters 
to the GUT parameters H IM E S H3 • Typically, this 
induces relations between the relative variations of the 
standard model parameters AGk/Gk- As an example, 
we estimate in Section^the variation AYne/^He for two 
models for the connection between the AGk/Gk within 
GUT models. The size of the effect depends strongly on 
these models. Keeping the ratio between the character- 
istic scales for the weak and strong interactions fixed we 
obtain AYjje/^He = 35.0 Aa em /a em whereas for a fixed 
ratio between the weak scale and the GUT scale one finds 
AFhc/1-Hc = 130 Aa em /a 

em • 



II. LINEARIZATION 

The success of BBN motivates the basic assumption of 
this paper, namely that the relative time variation of the 
fundamental constants between nucleosynthesis and the 
present epoch is small. We can then linearize in the rel- 
ative variation of the fundamental parameters AGk / Gk 
and use for Gk the values extracted from laboratory ex- 
periments. We express the relative change of the helium 



abundance as 

AYkc _ Y Bc (G + AG) - y Hc _ ^ (G)AGfc 

Ike >Hc ^ Cfc Gk ' 

Here Yho corresponds to the helium abundance computed 
in absence of a cosmological time variation of couplings, 
assuming that only standard model particles (with three 
neutrinos) contribute to the energy density at BBN. Our 

(G) 

aim is to determine the coefficients c k which relate the 
fundamental parameters to the change in Ye e - 

In this paper we will consider the effects of the varia- 
tion of six dimensionless quantities 

G k = ( M P /A QCD , 

&em i 

m e /A QCD , m q /A QCD , Am/A QCD ). (2) 

Here Aqcd is the characteristic mass scale of the strong 
interactions which dominates the mass of the nuclc- 
ons and the strong interaction rates whereas (4>) is the 
Fermi scale (vacuum expectation value of the Higgs field) 
relevant for the weak interactions. The strength of 
the gravitational interactions is given by the (reduced) 
Planck mass Mp and m e is the electron mass. The 
up- and down-quark masses m Ul are reflected in 
m q = (m u + m c i)/2 and Am = m ( i — m u . In combina- 
tion with (cj)}/AQCD the three last mass ratios could be 
replaced by the relevant Yukawa couplings h e ,h u ,hd- We 
emphasize that only ratios of mass scales are observable 
and have cosmological significance 0, IE H3 ■ 

For a given model for the time variation of the fun- 
damental constants the variations AGk/Gk are typically 
related to each other. For example, we may assume a uni- 
fied theory (GUT) and vary only the gauge coupling at 
the unification scale Mqut, while keeping (^3,4,5,6 fixed. 
This results in 5](AG 3 , 4 ,5,6 = 0) 

A{Mp/ Aqcd) _ tt Aa em _ tt a^ N - a em 
Mp/ Aqcd H a 2 em 11 a 2 em 

Then only a single independent varying coulping is left 
that we may choose as Aa em /a em . 

For practical reasons we will work in a frame where we 
keep the strong scale Aqcd fixed. This can be achieved 
by an appropriate Weyl scaling 0, 1^3] and will result 
in a time dependence of the reduced Planck mass Mp. 
This particular frame can be understood as a rescaling of 
the cosmological "clock" Mp which compensates for the 
constant strong interactions. In a frame with fixed Mp 
the strong interaction scale Aqcd would vary with time. 

As mentioned before our computation of the coeffi- 

(G) 

cients c\ ' proceeds in two steps. We first consider the 
dependence of Yhc on the characteristic quantities for 
nuclear decays and reactions, also referred to as "nuclear 
physics parameters" 

Xi = (Mp, a em , (0), m e , t„, Q, B d ), (4) 
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according to 



III. HELIUM ABUNDANCE 



AY, 



He 



He 



E 



AX, 

X, 



(5) 



Here, T n is the neutron lifetime, Q the neutron pro- 
ton mass difference and Bd the deuteron binding energy. 
(We keep Aqcd fixed - otherwise the dimensionful pa- 
rameters have to be multiplied by appropriate powers of 
Aqcd-) We emphasize that at this stage the effect of 
the variation of, say, a em is computed at fixed values of 

(X) 

^1,3.4,5,6,7- The computation of the coefficients c\ in- 
volves the details of nuclear physics, i.e. reaction rates 
etc. . 

A second step translates A Xi / Xi into the variation of 
the fundamental couplings 



2 j J ik Q 



Xi 



with 



fik 



d InXi 
d InGi. 



(6) 



This step involves the connection between nuclear physics 
and particle physics, namely the dependence of r„, Q 
and Bd on the couplings 6*2,3, 4, 5, 6- (Obviously, one has 
fik = $ik for i — 1...4 and k = 1...6.) For known fa- the 
coefficients c[ G ^ follow from X 



a.s 



(8) 



Before proceeding to estimates of the various coeffi- 
cients we list our results for the dependence of the helium 
abundance on the nuclear physics parameters (Section 
Unl in Table D 



TABLE I: Coefficients g for nuclear physics parameters 



variable 


Mp 


Bern {</)) rn e T n 


Q B D 


coeff. 


-0.81 


-0.043 2.4 0.024 0.24 


-1.8 0.53 



The dependence on the particle physics parameters is 
shown in Table [H] This Table allows for a quick inspec- 
tion of the impact of various parameter variations. We 
observe a particularly strong influence of a possible vari- 
ation of (0) and Am 



In this section we want to derive a semi-analytic esti- 
mate for the primordial helium abundance Yhc- In doing 
so we will follow the approach of Esmailzadeh, Starkman 
and Dimopoulos [38j], hereafter ESD. They estimate the 
BBN abundances via quasi static equilibrium and fixed 
point conditions. This approach should be sufficient for 
an understanding of the effect of small variations. Of 

(X) 

course, the estimate of the coefficients c\ as well as the 
determination of the corresponding errors would benefit 
from a systematic numerical investigation using the BBN 
codes. 

In the primordial universe, at temperatures above sev- 
eral MeV, the abundances of protons and neutrons are 
in thermal equilibrium. Protons are converted into neu- 
trons and vice versa. The latter process has a reaction 
rate given by |3^ | 



A \ dx x 2 f 1 



(Q + xf 



{Q + x) 2 , 

(l + e WT))-i(i + e -(Q+»=)/T)-i. (9) 

The integral runs from -co to +00 with an energy gap 
between —Q — m e and — Q + m e with Q being the proton 
neutron mass difference. Here A ~ {4>)~ i is the 4 point 
transition probability in Fermi-theory which depends on 
the axial and vector couplings cy and ca- For simplicity 
we will work with constant cy and ca- 

We now assume that this reaction freezes out at a tem- 
perature T* when the Hubble expansion is comparable to 
this reaction rate, i.e. 

T n ^ p {T*) = bH{T*). (10) 

Note that an equally well justified assumption would be 
to include the reaction rate T p -> n in this condition. The 
insufficient information of our simple approach (as com- 
pared to a more complete treatment by a solution of the 
Boltzmann equation) is accounted for by the unknown 
factor b. At the end we fix b so that we obtain the same 
Yhc as predicted by full numerical codes . The Hub- 
ble parameter H is given by the Friedman equation for a 
radiation dominated universe 2 



H 2 = 



P 



3Aff ' 



k "4 
9 * 30 T ' 



(11) 



with an effective number of degrees of freedom = 10.75 
before positron-electron annihilation. The freeze out 



TABLE II: Coefficients c k for fundamental couplings 



variable Mp a em (<j>) m e m q Am 
coeff. -0.81 1.94 3.36 0.389 -1.59 -5.358 



1 For b = 1 we obtain a 4 He abundance that deviates by about 10 
percent from the value Vjj e = 0.2484 found with a fully numerical 
computation using the WMAP value for r\ 12711 . In order for our 
analytic approximation to yield the Yfie predicted numerically 
we use b = 1.22. 

2 We do not consider changes in the expansion rate due to changes 
in baryon or electron masses. 
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temperature of the neutrons, with no change in funda- 
mental couplings, obtains as T* = 0.77 Mev and the 
neutron concentration at freeze out can then be calcu- 
lated as 



creating and 



Y * 



1 



,Q/T* 



0.158 



(12) 



Following the "freeze out" of the neutron to proton 
ratio the neutrons decay, thereby further changing Y n 
for T < T*. After a short time the synthesis of deu- 
terium and tritium starts which subsequently leads to 
the production of helium. Since almost all existing neu- 
trons end up in helium, we need to know how many neu- 
trons remained when helium was synthesized in appre- 
ciable amounts. We will assume that the neutrons decay 
freely until a time tf when helium formation starts to 
dominate over the neutron decay process, i.e. 

2Y Ue (t f ) = -Y n (t f ). (13) 

The final 4 He abundance is then estimated by 

Y He = \Y n {t f ) = Iy n * e -(*/A»>. (14) 

It depends on the couplings via Q, T*, r„ and tf. In 
turn, T* depends onA~ (c/))^ 4 , Q, m e and Mp via Eqs. 
©, flu and JTTJ)- 

We need an estimate of tf. The by far dominant pro- 
cess for helium production is the reaction ^(j 



d + t 



TIe + n 



(15) 



To write down the equations governing the abundances 
comprised of several reactions we will adopt the notation 
of ESD who abbreviate a reaction rate 



(16) 



as [a/37<5]. The condition for the time until which the 
neutrons decay is given by 



2Y d Y t [dtna] = —Y^f' 1 



(17) 



To compute the time when this relation is satisfied 
we need to know the abundance of deuterium (Yd) and 
tritium(Yt) as well as the reaction rate [dtna]. In the 
temperature range we are considering, deuterium can be 
assumed to be in thermal equilibrium and hence its abun- 
dance is given by the Saha equation 



Yd = 8.15 



T 



3/2 



77 e B *' T Y n Y p , (18) 



with the proton abundance being Y p « (1 — Y n ), m n the 
neutron mass. 

The estimate of Y t is more involved and also requires 
knowledge of the abundance Y3 for 3 He. The tritium 
concentration is established by the reactions 



s He+ n 
d + d 



p + t 
p + t 



(19) 



t + d 



l Rc + n 



(20) 



annihilating tritium. Other reactions are subdominant 
by at least 2 orders of magnitude (as can be verified 
from 01) an d are therefore neglected. Close to ther- 
mal equilibrium the fixed point condition leads us to 
an equation for Y t of the form: 



Y 



Y D Y 3 [n3pt] + Y d Y d [ddpt] 



(21) 



Likewise, we can write down the dominant processes for 
the 3 He abundance. Invoking the fixed point condition 
yields 



Y, 



Y d Y p [pd3 7 ] + Yd Y d [ddn3] 
Y3 [dipa] + Y n [nipt] 



(22) 



From Eqs. (JTHJ, $ZT\i and (|2*2|) we can determine the 
abundance of deuterium, tritium and helium-3 as a func- 
tion of T and Y n . In turn, temperature and time are re- 
lated by the background cosmology and Y n — Y* e~ l l Tn . 
Eq. (|13fl now determines tf. 

The dependence of Yh c on the various parameters can- 
not be solved analytically. In the linear approximation, 
however, the computation of the response coefficients 
c^' is straightforward. For this purpose we assume that 
all strong interaction rates are determined by the strong 
interaction scale Aqc d ■ At this point we benefit from our 
particular frame with constant Aqcd which implies that 
we can use constant rates [citna] etc., except for small 
electromagnetic effects. 

The results of this computation can be found in Ta- 
ble HJ They are plausible in the sense that they resemble 
what one would expect from simple arguments. Increas- 
ing the Planck mass gives a slower expansion rate, result- 
ing in a later freeze out of weak interactions, hence less 
neutrons are available for helium production. Increas- 
ing the decay time r„ of the free neutrons leaves more 
neutrons to be converted into helium since effectively all 
neutrons end up being bound in helium. Increasing {<j)) 
results in a decrease of the Fermi interaction Gf, hence 
weak interactions freeze out earlier resulting in an in- 
crease in Yho- Changing Q results in a different neutron- 
proton ratio at freeze out and also in modified weak rates 
due to changes in the available phase space. If we exclude 
the changes in phase space volume, the coefficient is —1.4 
instead of —1.8. Thus, helium abundance is a decreasing 
function of the proton-neutron mass difference Q as an- 
ticipated. Increasing the binding energy of the deuteron, 
Bo, results in earlier formation of helium and reduces 
the amount of neutrons decaying into protons. The influ- 
ence of the electron mass is only through the phase space 
volume in the weak rates which is a very small effect for 
our purposes. 

Changes in a em affect the nuclear reaction rates with 
the main effects being variations in the Coulomb bar- 
rier for charge-induced reactions, final-state interactions, 
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radiative capture and mass differences. We use t he p ro- 
cedure of Bergstrom, Iguri and Rubinstein (BIR) [43] for 
computing the impact of varying a em on all nuclear re- 
action rates used in our computation, including the im- 
provements of Nollett and Lopez |4^ | . We use the rates 
of the NACRE compilation [44J where available, other- 
wise we use those of Smith, Kawano and Malaney (SKM) 
[icj . For the process 3 He(n,p)i we use the fit of Cyburt, 
Fields and Olive . Since the analytic NACRE rate fits 
have a different expansion in terms of T we have fitted 
the rates to the SKM functional form for use of the BIR 
treatment as described in |43|. 

Except for electromagnetic effects we have not taken 
into account any other effect that may change the reac- 
tion rates. 



IV. FUNDAMENTAL COUPLINGS 

In this section we describe the relation between the 
fundamental couplings Gk and the nuclear physics pa- 
rameters Xi. This relation was expressed in the form 
of a matrix equation J7J). We will now discuss what ef- 
fects we took into account by explaining each row of the 
matrix fa (see Tabic ILTlj) . Each coefficient fik describes 
the response of the "nuclear physics parameter" Xi when 
one varies a single parameter Gk , while keeping the other 
Gi-tk fixed. For i — 1...4 the parameters appear both in 
the lists of Xi and Gk and fa = Sik by virtue of our 
definition. Also r n , Q and Bd do not depend on Mp 
implying fik = 5\k- The nontrivial coefficients fik for 
i = 5,6,7 account for the dependence of r n , Q and Bd 
on a em , (<f>), m e , m q and Am. 

The nucleon masses and nuclear binding energies de- 
pend on the quark masses and a em . The dependence of 
the neutron-proton mass difference on the fundamental 
couplings is given by (see pfij): 



Q = 



-0.76 1 



Ac 



2.05 1 



A(Am) 
Am 



MeV . 

(23) 

From this we can determine fe2 and fee- Recent studies 
have suggested that the deuteron binding energy Bd may 
increase with decreasing pion mass 0, E3- We may 
parametrize the dependence of Bd on m^ at fixed (<fi) by 
a linear fit |r| and neglect the dependence on (0) at fixed 
m q , Am. For the electromagnetic part we use the Monte 
Carlo simulation data of Pudliner et al. [Sjj. Hence the 
deuteron binding energy may be expressed in terms of 
the pion mass m n oc vn'J' 1 and a em as 



Bd - B d 



(r + 1) -r— 

rrvi 



0.018^^ MeV, (24) 



where r is a parameter that varies between 6 and 10 
and B® — 2.225 MeV is the deuteron binding energy as 
measured in the laboratory today. 

The neutron lifetime is changed due to variations in 
the weak scale r„ tx Gp 2 oc ((f>) 4 . Furthermore, a change 



in the phase space volume / of free neutron decay 
rQ 



f = 



dq q'(Q-qY(l- , 
Q 



e\l/2 



(25) 



results in a dependence of r„ on Q and the electron mass 
m e . Because Q also depends on a em and Am, Ar„ will 
also have contributions from the variation of those pa- 
rameters. A linear analysis then yields the corresponding 
entries for r„: 



Ar„ 



= 3.86- 



Aa. 



A(d>) Am e A Am 
4— ^p- + 1.52 - 10.4—^ '-. 

T n a em (<p) m e Am 

(26) 

The entries of the matrix fa can be found in Table 
IIIII We have quoted the coefficients for the effects we 
discussed above. When there is no contribution at all 
a zero is written. For some relations between the Gk 
and the Xi small effects are present but with negligible 
coefficients. To distinguish those from the others, we 
have left the matrix entry empty. Having determined 
the transfer matrix we can calculate the dependence of 
AYne/^He on the fundamental parameters (see Eq. JSJl). 
The results are shown in Table iHl above. 



TABLE III: The matrix entries fik, corresponding to the 
coefficients relating relative changes in Gk to relative changes 
in X t . 



parameter a 


Mp 




w 


m e 


m q 


Am 


Mp 


1 






















l 






















1 






















1 








T„ 





3.86 


4 


1.52 




-10.4 


Q 





-0.59 








1.59 


B d 





-0.0081 






-r/2 





"The parameters are dimensionless, but we omitted the scaling by 

QCD 



A 



V. TWO GUT EXAMPLES 

In this section we present two examples based on GUTs 
for which we have expressed the changes in the funda- 
mental parameters by the variation of only one indepen- 
dent coupling and computed the resulting change in Ihc- 
The variation of the couplings is assumed to be due to a 
scalar field y called the cosmon There are good 

arguments p that this field plays the role of quintessence 
|30j | and its present potential and kinetic energy can be 
associated with the dark energy of the universe. For our 
considerations, however, we will not need any particular 
details of the evolution of the cosmon, except that its 
value at the time of nucleosynthesis was different from 
its present value. 
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For the details of the derivation of how the fundamen- 
tal constants change in a GUT scheme we refer the reader 
to Merely quoting the results, to one loop order 

the fundamental couplings as functions of the cosmon 
field x are given by [5| 



— - 2 — + Tr^Cw{x), 

9 2tt 
^Z F ( X ) , 5 



r 



3?r 



m C™(x), 



32?rZ F (x) 5 

— gps— "a^to, 



(27) 
(28) 
(29) 



where the W-Boson mass is Mw(x) — Cw(x) X and 
Zp(x) determines the renormalized grand unified gauge 
coupling (g 2 R = g 2 /Zp, g fixed). We normalize x such 
that Mcut(x) = X- In Eqs. $Z7$-&$ we can replace 
Mw — 9to(0)/v2 by {(f)). The relative variation of g w 
(or a w ) induces only a correction of higher order in these 
relations. 

As mentioned before, we will work in a frame in which 
the scale of the strong interaction is fixed such that the 
strong interaction rates are constant for our BBN esti- 
mate. We will consider two particularly simple scenarios 
where Mp(x)/Mgut(x) — const, with 



AMp/Aqcd 
Mp/AqcD 



-Alnfu, + Aln.({<f>) / Aqcd)- (30) 



Furthermore, we also neglect the variation of the 
Yukawa couplings and hence the variations in m e , m q 
and Am obey 



A? 



A(Ato) 



Am„ 



(0> ' 



(31) 



The effect of the variation of the cosmon field x can now 
be expressed as a varation in the renormalized grand uni- 
fied gauge coupling expressed by Zp and a variation in 
lnCtu- 

At this stage the two unknown quantities A In Zp and 
A In contain all relevant information about the un- 
known coupling of the cosmon to matter and radiation. 
For the present investigation we can simply use the rela- 
tive variation of the GUT-coupling A In Zp and the ratio 
between weak and GUT scale A In ( w as free parameters. 
For the running of a em at /i < Mw we have the relation 



,(Af! 



w) 



3- 



In 



M w 



(32) 



where the Qi are the charges of the particles with masses 
in the range between My/ and fi. In our case this is given 
by five quarks (top lies above My/) in three colours plus 
3 leptons, i.e. £\ Q? = 3 x (8/9 + 3/9) + 3. Similarly, 
for the running of a s below My/ we include five quarks 
and associate Aqcd with the scale where the one loop 
expression for a s (p) diverges. 



We can now express a em = a em (m e ) and Aqcd in 
terms of a em {My/) and a s (My/). Thus they relate 
Aqcd/x an d a em to Zp and lnCw- The specific relation 
between the variations of Aqcd/M p and a em depends on 
the variation of the weak scale In ( w . Our two examples 
will either keep ln^ or ln(((f>) j Aqcd) fixed. 

The first example is as simple as possible - we also 
keep the weak scale fixed w.r.t. the strong scale. This 
will result in a x independent ratio {(f)) / 'Aqcd ~ const, 
and leads to 



Aa em {M ] 



w) 



in AZ F 



and 



AMp/Aqcd 

Mp/AgcD 



7 



7r Aa e 

TT 



.{My/) 



(\ - 



,(M, 



(33) 



(34) 



Since the weak scale is fixed we set A{4>)/{4>) = in 
Eq.J]J. 

In our first example one has Aa~^(m e ) = Aa~^{My/) 
which is related to the only unknown parameter AZp by 
Eq. ■ Eq. $Q results in 



A{Mp/A QCD ) = _ 3g 1 Aa em 
Mp/Aqcd 



(35) 



The only nonvanishing entries in Eq. from AG± 

and AG2 are therefore related and AYn e /Yn e can be ex- 
pressed in terms of a single parameter that we may choose 
as Aa em /a ern . In order to get an idea of the sensitivity 
we compute the value Aa em /a em which would be needed 
in order to obtain a helium abundance Yho = 0.24 for ij 
corresponding to the central WMAP value. For our first 
example we find 3 



Aa em (m e ) 
a em {m e ) 



-1.0 x 10" 



(36) 



In the second example we will set {(f>) jx — const. For 
this setting we obtain 



Aa em (M w ) 
a 2 em {M w ) 



32tt AZ F 



and 



AM p/ Aqcd _ A{(f>) it Aa em {M w ) 



M P / Aqcd 



(0) 



12 al m {M w ) 



(37) 



(38) 



We also need to relate the different mass scales via 
Eq.(j23) giving 



Ac 



(m e ) _ Aa 

cm- 

(M w ). 



(M w ) 



(39) 



3 We have put r = 6. The change in Yg_ e is negligible if we choose 
r = 10. 
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where J^iQi = 2 runs only over the three light quarks 
whose effect on the running of a em is cut off at fi ~ 
Aqcd- 

In Eq.Q) we now have nonvanishing entries from 
^^3,4,5, 6 as well, related by Eqs. (joT]) and to 
Aa em /«em giving 



A(M p /A QC d) _ A(<j>) 



Mp/Aqcd 



<</>> 



-32.3 



Aa e 



(40) 



The variation of Aa em which would be needed for 
Yke = 0.24 is now reduced as compared to the first model, 



Aa -} m ^ = - 2 .7 x W-\ 
a e m{rn e ) 



VI. DISCUSSION AND CONCLUSION 



(41) 



Big Bang nucleosynthesis offers an excellent testing 
ground for the time variation of the fundamental cou- 
plings, probing directly their values at a time close to 
the big bang or at high redshift. This paper presents a 
general analysis how the primordial helium abundance 
depends on six "fundamental couplings" as summarized 
in Table [H] In grand unified (GUT) models the varia- 
tions of the various couplings are interrelated. Typically, 
the dominant effect comes from a variation of the ratios 
Mp /Aqc d or (0) / Kqcd rather than from the direct in- 
fluence of a variation of the fine structure constant p|- 
This is easily seen by comparing Eqs. (|35[1 and H4U|) with 
Table [n] As noted before |5| the variation of the weak 



interaction scale (<f>) or the quark masses can play a very 
substantial role. We have checked that the impact of a 
variation of the strong and electromagnetic interaction 
rates is a rather minor effect. Discussing two different 
models gives us a certain handle to investigate the effect 
of cancellations - compare Eqs. (|3u[l and l|41|) . Extending 
our analysis to the abundances of deuterium and lithium 
may be a way to (partially) lift the degeneracies between 
the variations of various couplings. 

Excluding very particular cancellations we may in- 
fer from the approximate agreement between the 
WMAP-prediction and the observations of Yho a bound 
\Aa em /a em (z = 10 10 )| < a few times 10~ 3 . A typical 
size of a coupling variation that could explain the present 
discrepancy between WMAP and the observed Yhc would 
be in a range Aa em /a em » (2 — 10) x 10~ 4 . 

Obviously our treatment can be improved. The coeffi- 
cients quoted in Tabled could be estimated with higher 
accuracy by using a full numerical code instead of our 
analytic estimate. Likewise, looking at Tabic ITTT1 one can 
see that there are some small effects contributing to the 
matrix fa- which we have not included. Also, we have 
only investigated the change in the helium abundance. 
Stringent bounds on BBN are obtained from the primor- 
dial abundances of deuterium and lithium and it would 
be worth extending our analysis to these other light ele- 
ments. 
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